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Abstract. The isomonodromy deformation method is applied to 
the scaling limits in the linear NxN matrix equations with rational 
coefficients to obtain the deformation equations for the algebraic 
curves which describe the local behavior of the reduced versions for 
the relevant isomonodromy deformation equations. The approach 
is illustrated by the study of the algebraic curve associated to the n- 
large asymptotics in the sequence of the bi-orthogonal polynomials 
with cubic potentials. 



1. Introduction 

It is well known that, in certain asymptotic limits, the classical 
Painleve equations [1] reduce to elliptic ones. For instance, the Painleve 
sixth equation, 

1/1 1 1 \ 2 (\ 1 1 \ 

yxx = 7;(- + — 7 + )yl-{- + — 7 + )Vx+ 

2 \y y — 1 y — x/ \x x — 1 y — xJ 
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with the large parameters Cj, j = 1, . . . , 4, after the changes x = to+5r, 
Cj = 5~ 2 cij + $~ l bj ) where t Q ,aj,bj = const, turns at 5 = into an 
autonomous equation which has the first integral 

D o = 7Tl 777 -ry T - a iy + a 2- + a3 7- + °4 — . (1.1) 

2y(y- i)(y-*o) y y~ l y-h 

The asymptotics of the classical Painleve transcendents w.r.t. parame- 
ters or initial data were studied in numerous works, see [2] for extended 
but not exhaustive bibliography. The limit transitions of such kind are 
called below the scaling limits. 

The most effective to this date approach to the scaling limits in the 
Painleve transcendents is based on the monodromy representation for 
the latter [3, 4]. In [2], the isomonodromy deformation technique of [5] 
was adapted to the study of the scaling limits in the equations of the 
isomonodromy deformations for the 2x2 matrix linear first order ODEs 
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with rational coefficients. In particular, the results of [2] imply that 
the modular parameters determining the limiting (hyper)elliptic curve 
for the asymptotic solution of the isomonodromy deformation equation 
like D in (1.1) are not arbitrary constants but certain functions of all 
the deformation parameters. These functions are uniquely determined 
by the system of transcendent modulation equations 



where £ is an arbitrary closed path on the Riemann surface of the 
relevant spectral curve. 

Below, this result will be extended to the scaling limits in the isomon- 
odromy deformation systems forNxN matrix linear ODEs with ratio- 
nal coefficients. The work is motivated by recent developments in the 
theory of coupled random matrices. Indeed, while the statistic prop- 
erties of the ensembles of the single random matrices can be given in 
terms of the asymptotics of the semi-classical orthogonal polynomials, 
see [6, 7], which give rise to the linear first order 2x2 matrix ODEs 
with rational coefficients [8] , the ensembles of the coupled random ma- 
trices in the very similar way give rise to the bi-orthogonal polynomials 
[9, 10, 11] and the linear N x N matrix ODEs [10]. 

The paper is organized as follows. In Section 2, we recall the basic 
facts in the isomonodromy deformations of the linear matrix equations 
with rational coefficients, introduce the notion of the scaling limits 
in such systems and describe the WKB approach to their asymptotic 
solutions. In Section 3, we find the modulation equations for the non- 
singular asymptotic spectral curve and prove their unique solvability. 
In Section 4, we illustrate our approach using a particular 3x3 ma- 
trix equation satisfied by the bi-orthogonal polynomials for the cubic 
potentials. 

2. The isomonodromy systems with a large parameter 

In this section, following [12, 13, 3], we recall the basic notions of the 
theory of the linear matrix differential equations. Consider an N x N 
matrix first order ODE, 



JVTi 71 Tv . Too 

f x ^A( m aw = ee ^_- Ea ^-k (2 ,) 

v=\ k=0 v ' k=l 



We call equation (2.1) generic if the eigenvalues of A v _ Tv are distinct 
for r„ ^ and if they are distinct modulo integers for r v = 0. Without 
loss of generality, A = oo is the singular point of the highest Poincare 
rank, i.e. > r v , v = 1, . . . , n. Assume that (2.1) is generic and 




(1.2) 
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Aoo- roc is diagonal (for the non-generic situations, see [13]). Then, 
near the singularity aS v \ equation (2.1) has the formal solution 

tf„(A) = W {v) ¥ v \\)e e{v){x \ v= l,...,n,oo, 

3=1 3=1 ' J ' 

where £ = A — for a finite singularity and £ = 1/A for infinity. 
The matrix coefficients ip^ and the diagonal matrix coefficients x^j of 
the formal expansion (2.2) are determined uniquely by the eigenvector 
matrix of A v - Tv . 

The ratio ^ _1 (A)^(A) of any two solutions ^ and \& of (2.1) does 
not depend on A. The ratios of the fundamental solutions normalized 
by (2.2) are called the monodromy data. 

The set of deformation parameters is specified for generic equation 
(2.1) in [3] (for non-generic equations in [14]). These parameters in- 
clude the positions of the singular points and the entries of the 
diagonal matrices x^j, j ^ 0, for ^(A) in (2.2). All these quantities 
form together the vector x of the deformation parameters. Remain- 
ing parameters x^\ v = 1, . . . , oo, are called the formal monodromy 
exponents. The latter satisfy the Fuchs' identity, 

n 

^Tr^+Tr^ = 0. (2.3) 



Remark 2.1. Using the linear transformations of the complex A-plane, 
one can fix the positions of two of the finite singular points or to fix 
two of the entries (x^) kk . Using the scalar gauge transformation, one 

can get Tr A(X) = and therefore Trx^j = 0. The orbits of the above 
transformations are called the essential deformation parameters. 

Let d denote the exterior differentiation w.r.t. entries of x. In accord 
with [3], the monodromy data of the generic equation (2.1) do not 
depend on x if and only if there exist 1-forms Vt and 6*^ such that the 
fundamental solutions above additionally satisfy equations 

dV = m, dW {u) = <d {u) W {u \ u = l,...,n,oo. (2.4) 

Then the compatibility condition of (2.1), (2.4), 

dA = —+[Q,A], dQ = QAQ, (2.5) 
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is the completely integrable differential system whose fixed singularities 
are the planes = a^ p \ v ^ p, = oo, {x^} rv )u = 0, {x^} rv 



in 



(x-ljjj if * 7^ j, r u 7^ 0. The generic 2x2 system (2.5) admitting the 
only one deformation parameter is equivalent to one of the classical 
Painleve equations. 

2.1. Scaling limits in the isomonodromy systems. Let A(X) de- 
pend on an additional parameter 8, 

a H = S x, b (y) + S (u) \ const = 1 

(2.6) 

A v -k = 5 K B v _^, ^4oo,-fc = ^ x ^oo,-fe) /c = 0, . . . , r„, 

where B v _ k , k = 0,...,r u , v = l,...,n, oo, are ascending Erdelyi 
series in <5 such that, for generic equation (2.1), 

o>) = + St™), = + 5r5 } ). (2.7) 

The constants and the entries of t^\, k ^ 0, form the vector £ for the 
center of the asymptotic domain in the parameter space as 5 — »• +0. 
Inequalities |c^|, ||t^|| < const, or simply ||r|| < const, yield the 
range of the "local" deformations. The entries of t are usually called 
the "slow" variables, while the entries of r are called the "quick" or 
"fast" variables. The entries of the vectors t^ and for the formal 
monodromy exponents form the vectors a and f3, respectively. 
Substituting (2.6) and A = 5"( into (2.1), (2.4), we find 

^ = 5- 1 B(C)^, d* = u*, (2.8) 

Tl Ti/ j-^ Too 

v=l k=0 VS ; k=l 

whose compatibility reads 

dB =[uj,B} + 5^, du = uf\u. (2.9) 

Remark 2.2. If ^ and (x^Jn does not depend on 5 then x = 
— I/Voq. If all the singularities are Fuchsian, i.e. r„ = 0, v = 1, . . . , n, oo, 
one may put for simplicity h—0. 



Remark 2.3. Following [15, 8], it is possible to construct a Schlesinger 
transformation chain yielding x^ = + with constant G 
and G Z*. Taking and 5" 1 large and comparable, 
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lim^+o^ 5 = = const, one writes = 5 + where 
Tq^ G C n is bounded as 5 — > +0, and 

t ( u) G R N . (2.10) 

This allows us to interpret Re {t^)jj as the discrete deformation pa- 
rameters with the step 0(S), see below. 

Remark 2.4. The scaling parameter 5 is defined up to a positive factor 
which gives rise to a scaling freedom in the set of the "slow" variables. 
Below, we assume that this scaling freedom is eliminated by a normal- 
ization of one of the non-trivial "slow" deformation parameters. 

2.2. Complex WKB method. Here we recall the idea of the complex 
WKB method following in principal [12, 13]. Consider (2.8) as 5 — > 
assuming that the coefficients of B(() remain bounded. Let T and A 
be the eigenvector and eigenvalue matrices for B((), 

T- 1 J BT = A = diag(/i 1 ,..., / u JV ). (2.11) 

Let \ij 7^ /ifc, j 7^ k. Then the formal expression 

OO OO 

y(() = Tj2t n T n exp{5- 1 / ^5 m A m (0^}, T = I, (2.12) 

ra=0 m=0 

satisfies (2.8) provided the diagonal matrices A n and the off-diagonal 
matrices T n , n > 1, solve the recursion 

[A^J-A^f, 

[A„, r„] - A„ = ]P X„.„,A m + T' l d{T f'- l} ■ n>2. (2.13) 

m=l 

Let T be the Riemann surface of the algebraic curve 

F(C,/i):=det(5(C)-^/)=0. (2.14) 

The branch points of (2.14) are called the turning points for (2.8). Let 
Cq be an open simply connected domain which is the complex (-plane 
punctured at the singularities of B((), at the turning points and cut 
along the segments connecting all the singularities and turning points. 
Let C C Co be a closed simply connected domain. By construction, 
B(() is holomorphic in C, and equation (2.8) has no turning point in 
C. Thus all the roots fJ>j(X) °f the characteristic equation (2.14) are 
distinct from each other and are holomorphic in C. Therefore there 
exists a holomorphic non-special in C matrix T(Q which diagonalizes 
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the matrix B(Q, T~ X BT = A = diag(//i, . . . ,/j,n)- The matrices T, 
T _1 , A n , T n , are holomorphic and bounded in C [12]. 
Consider the reduced gauge matrix 

m 

T^=Tj2$ n Tn, T = I, (2.15) 

n=0 

where T n are defined by (2.11)— (2.13), and ( G £. The matrix function 
$(m) = (7-(m))~i^ solves the "almost diagonal" equation 

m 

$( m ) = S^BWQW, B {m) {() = ^<TA n + 5 m+1 J R (m) , (2.16) 

?1=0 

where R( m \() is holomorphic and bounded for £ G C provided 5 is 
small enough. Define the WKB approximation to (2.16), 

K m 

"'Co n=0 

(2.17) 

and introduce the correction function x^ m \ 

$(-)(C)=x (m) (C)^£ B (C), (2.18) 

and notations 

^ = ^ - %(e, o = we, o)« - we, o)„., 

"C (2.19) 



The contour 7i 3 -(C) C C connecting the finite or infinite point Qj with 
( is called the (i, j)- canonical path if 



Re J C finals <0 V£ G 7i?(C)- 



(2.20) 



A closed simply connected domain C C is called (i, j)- canonical 
if there exists such a point ^ G that the contour 7ij(C) C 
connecting £y with any given point ( G Cij is homotopy equivalent to a 
canonical path. A closed simply connected domain C is called canonical 
if it is (i, j)-canonical Vi, j = 1, . . . , N, i ^ j. 

Using the arguments of [12], we thus obtain the following 

Theorem 2.1. Let C C C be a canonical domain. If 
[ e Rc§ ^\{R^(0)J-m<oo, Vi,j = l,...,N, VCGC, 
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then there exist such positive constants C and 5q that 

\\x {m ' 1] (0- I\\<C\5\ m VCeC, V5e(0,5 ]. (2.21) 
In particular, x (0) (C) = I + 0(6). 

To construct a canonical domain containing a given point Co £ £, 
consider a pair of i, j G {1, . . . , iV}, % ^ j, and introduce the segment 
£ij C £ of the (i, j) -anti- Stokes level curve-line passing through ( : 

£i j = {(e£: Im / ^-(s)ds = 0}. (2.22) 
Jco 

Choose two points QjXji £ m suc h a way that: a) £ £ [Cr>?Cji]; 
b) for any £ G % separating ^ from Q iy the curve-line segment 
[Cij,C] C % is the (i, j)-canonical path, while [Cji>C] C % is the (j, i)- 
canonical path. Introduce the segment £*j C £ of the (i, j)-Stokes level 
curve-line passing through the point (* G [Cij? 0*] c 

£*. = {( e £: Re / C ^(s) ds = 0}. (2.23) 

By construction, the union dj of all the curve-line segments £*j, 

Ca= U ( 2 - 24 ) 

is the and (j, i)-canonical domain. The boundary of the con- 

structed (i, j)-canonical domain dj is formed by the (i, j)-Stokes level 
curve-lines passing through the points Qj and Qi and partially by the 
boundary of £. It is worth to note that, near the irregular singularities 
of (2.8), the (i, j)-canonical domain dj can be extended beyond the 
boundary of £ to fill out certain sector in the complex £-plane called 
the (i, j)-Stokes sector. 

The canonical domain C 3 Co of validity of Theorem 2.1 is the inter- 
section of the above (i, j)-canonical domains dj, 

C= p| dy (2.25) 

i,je{i,...,N} 

The above construction implies that any closed simply connected 
domain £ C £o can be covered by a finite number of the overlapping 
canonical domains d k \ k = 1,...,K, since the opposite assumption 
can be easily brought to a contradiction. 
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3. Modulation of the spectral curve 



Solutions of the Lax equation dB = [u, B] are routinely interpreted 
as the approximate solutions for (2.9) as 5 — > [16, 17]. Supplemented 
by the eigenvalue problem (2.11), the Lax equation constitutes the basis 
for the algebro-geometric integration of the "soliton" equations. How- 
ever, in the theory of the "soliton" PDEs, the spectral curve (2.14) is 
determined by the initial data, while in the isomonodromy deformation 
context it is determined by the original A-equation (2.1). Moreover the 
spectral curve for the typical "soliton" equation is an exact integral of 
motion while, in the isomonodromy case, the curve varies, 



In what follows, we precisely describe the dependence of the algebraic 
curve (2.14) on the "slow" variables at S — 0. Below, the subscript as 
denotes the relevant object at 5 = 0. 

We call the curve (2.14) singular if its topological properties for all 
small enough 5^0 differ from those at 5 = 0. Given a parameterization 
of the curve, we define the discriminant set & in the total parameter 
space ^3 as the set determining the singular curve. Also, let $ be the 
union of the hyper-planes b {u) = b ip \ v ^ p, (t^lju = (i^Jjj, % ^ j, 

(t-l v )kk — corresponding to the fixed singularities of (2.9) at 5 = 0. 
Below, we always assume that our deformation parameters are apart 
from the fixed singularities 

The differential fi as (() d( as well as its derivatives w.r.t. and 
entries of t^l, k — 0, . . . , r v , v = 1, . . . , n, oo, are meromorphic on the 
Riemann surface r as of the curve. All the parameters b^ , &\ together 
completely determine the singular part of fi as (C) d(- 

Definition 3.1. The parameter Dj is called modular iff the differential 
a§~A t as(C) is holomorphic on the Riemann surface T as . 

Thus = 1 <g> 2), where 1 is the subspace of the deformation pa- 
rameters t = (t,Rea) (see Remark 2.3 and constraint (2.3)) while D 
is the subspace of the remaining parameters D = (D, Ima). 

Theorem 3.1. Let (t ,D ) G ^P\<5. Then there exists an open neigh- 
borhood U C X\# of t such that, for any closed path t on the Riemann 
surface Y as punctured over the points b^ u \ v — 1, . . . , n, oo, 



dflndet (B -fjil)) = 5Ty(—(B - fiiy 1 ) ^ 0. 




Vt e U, 



(3.1) 



where he = J^(t , D ) = const. 
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Proof. Let £^ = n(£) be a projection of the closed path £ C T with the 
base point (Co,A*o) 011 the punctured complex £-plane. Let the integer 
m be chosen in such a way that the lift £ of m £^ on T be closed for 
all branches of //(Co)- Consider the analytic continuation of the WKB 
approximation (2.12) along £. The projection n(£) = is covered by 
a finite number of the overlapping canonical domains Ck, k — 1, . . . , s, 
C s+ i = Ci, in each of which (2.12) approximates uniformly in ( an 
exact solution ^k(C) of (2.8). Because — ^k(C)Gk where G k 

is independent from both ( and t G X \ we obtain 

M l (* 8+1 (0) = * 1 (QMp 1 ...G 8 , 

where M.$ is the operator of analytic continuation along mo^, and M| 
is the monodromy matrix for ^(C) along m £^. Using for \&i(C) and 
^s+i(C) our WKB approximation, we find 

expjr 1 jf A(C) d() = (I + (3.2) 

Since the curve is non-singular, the r.h.s. of (3.2) preserves while t 
remains in a neighborhood of to- Equating the leading orders of the 
l.h.s. for (3.2) at t and nearby points t, we arrive at (3.1). □ 

Theorem 3.1 immediately provides us with the following assertion: 

Corollary 3.2. Let U C X \ $ be an open domain and let (3.1) holds 
true. If the curve F (£,[/,) = remains non-singular at the boundary 
point ti G dU, then there exists an open domain Wcl\jf such that 
U CW, ti eW, and (3.1) is valid Vt G W. 

For the subsequent discussion, the following assertion is useful: 

Proposition 3.3. For any closed path £ on Y as punctured over b^ u \ 
v = l,...,n, oo 7 the integral J^(t,D) is continuous in (t,D) outside 
the fixed singularities $ of (2.9) and is differentiate in (t,D) outside 
the discriminant set &. 

Proof. Since the contour £ is finite, the continuity of J^(t,D) follows 
from the continuity of fi as ((). If the point (t ,D ) is located apart 
from the discriminant set &, then there exists an open neighborhood 
V of (t , D ) such that the spectral curve (2.14) does not degenerate 
V(t,D) G V. Then the differentiability of J*(t,D) at (t ,D ) follows 
from the continuous differentiability of /i as (C). □ 

Theorem 3.1 and Proposition 3.3 imply 
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Corollary 3.4. Let U,U C X \ $ be adjacent open domains and let 
(dU H dU) C 1\$ be not empty. If J e (t, D) = hp Vt G U and 
Je(t, D) = hp Vt G U, then he = hp. 

In accord with Corollaries 3.2 and 3.4, if the spectral curve (2.14) is 
non-singular at the initial point (t , D ), then the modulation equation 
(3.1) is valid in a domain U bounded by the points where the spectral 
curve becomes singular. Applicability of (3.1) to a particular solution 
of (2.9) beyond this boundary depends on some subtle details in the 
initial data or, equivalently, in the relevant monodromy data of the 
isomonodromy system (2.8), see [2] and Section 4 below. 

Let us discuss now the existence of the function D(t) such that 
j£(t,D(t)) = const. Varying the contour £ in (3.1), we obtain the 
system of equations = hj, where the set of contours £j form a 
homology basis of the Riemann surface of T as punctured over . For 
instance, taking for £ a small circle c v around (b^\ fj,j(b^)), we find 

Im(#)).. = -^g) = const, (3.3) 

which contains (2.10) as the particular case he] = 0. 

To discuss (3.1) further, it is convenient to impose the conditions 
(3.3) and to remove (N — l)(n+l) small circles from the "sufficient" set 
of contours. The remaining cycles £j, j = 1, . . . , 2g, form a homology 
basis of T as . Also, using (3.3), we exclude the constant parameters 
Imto^ from the set of unknowns and assume below that the space D 
is g-dimensional complex space of the modular parameters. 

Theorem 3.5. Let (t ,D ) G ^3 \ &. Then, in an open neighborhood 
U C X \ 5 of the point t 0; the system (3.1), where £ runs over the 
homology basis {£j}\ 9 of T as , determines the unique differentiable in 
the real sense complex vector function D{t, t) such that D(t , t ) = D . 

Proof. Theorem 3.1 and Proposition 3.3 imply that J^.(t,D) are the 
first integrals of the completely integrable Pfafnan system dJ = 0, 

»(S)~°(2). <-> 

where u and Q are the matrices of the partial derivatives of Jg {t,D) 
w.r.t. the entries of the vectors D, D and t,t, respectively. Here, the 
bar means the complex conjugation. Let constant C\ be the determinant 
of the transformation of the natural basis {jjj-[i as c?C}f=i into the basis 

of the normalized holomorphic differentials, and let B be the matrix of 
the .B-periods of the normalized holomorphic differentials. Since uj = 
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( g 4 ) is the matrix of A- and 5-periods of the holomorphic differentials 

and their complex conjugate, detu; = (— 2i) g \c 1 \ 2 det(lmS) ^ 0. Thus 
the matrix uj is invertible until F (£,[/,) = remains non-singular, and 
therefore the integral manifold for (3.4) is well parameterized by the 
deformation parameters (t,t, Re a). □ 

Remark 3.1. If = J^(t , -Do) 7^ then the cycle £j can not collapse, 
and the encircled by £j branch points can not coalesce. Thus, along 
the integral manifold for (3.4), the spectral curve remains non-singular 
provided ^ Vj = 1, . . . , 2g. This observation ensures the appli- 
cability of (3.1), (3.4) and the existence of D(t,t) in any connected 
domain WcI\S containing the initial point t . 

Let j£(t,D) in (3.1) vanish for all closed paths, 

Re jf// afl (C)dC = 0, WcL. (3.5) 

This system does make sense regardless the choice of the initial point 
since there is no need to fix a homology basis. Traditionally, it is called 
the Boutroux system. We recall that (3.5) may be not applicable to a 
particular solution of (2.9) in certain sectors of ^3 \ & in spite of the 
Boutroux system itself does make sense in the whole parameter space 
*P (using Proposition 3.3, the system (3.5) is interpreted at the points 
of the discriminant set © as a continuation from *}3 \ &). 

Remark 3.2. As the integral manifold for (3.5) meets the discriminant 
set (5, at least one of the cycles £j collapses, and the corresponding 
real equation in (3.5) becomes trivial as being replaced by the complex 
condition of coalescence of two branch points. Thus, generically, the 
intersection of the integral manifold for (3.5) with & has codim^ = 1 
in the space of the deformation parameters. 

Theorem 3.6. There exists the unique solution _D(t, t) of the Boutroux 
system (3.5). 

Proof. Here, we give the sketch proof. 

Uniqueness. Given t, two solutions D and D' determine two dif- 
ferentials jj, as d( and n' as d( meromorphic on the respective Riemann 
surfaces r as and T' as . The difference <fi = (fi as — fi' as ) d( is holomorphic 
on the covering Riemann surface Q as , and Re § t <f> — for all closed 
paths £ C Gas- However, there is no differential with such properties. 

Existence. Choose a point (t , D ) G *p\© in such a way that hf = 
J^.(t , -D ) 7^ for all contours £j of a homology basis {£j}\ 9 ■ Consider 
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the extension of (3.4) where 2g real parameters h — (hi, ... , h2 g ) T are 
added to the set of the independent variables, 



Applying the arguments used in the proof of Theorem 3.5 and taking 
into account Remark 3.1, we establish the existence of the function 
D(t,t,h) Vt e U C T\ 3 and V/i: /i i sgn(^ 0) ) > 0, j = 1, . . . , 2g. 
The assumption that D(t,t, h) is unbounded as h — > leads to a 
contradiction. From a bounded sequence = -D(t,t, h^), — > 0, 
we extract a convergent subsequence, linim^oo D^ km ^ = D*. Then the 
continuity of the integrals j£ 3 (t,D) w.r.t. D yields J^(t,.D*) =0. □ 

Remark 3.3. Assuming that the monodromy data of the system (2.8) 
are generic and do not depend on the scaling parameter 8 -1 , or this 
dependence is weak enough, it is possible to prove that the relevant 
spectral curve satisfies the Boutroux system (3.5). Here, we do not 
prove this assertion (look for more details and for the proof of this 
statement in the case N — 2 in [2]). 

4. Modulation equations and the asymptotics of the 



Modulation equations (3.1), (3.5) for the classical Painleve equations 
were studied in [18, 2}. Here we note that the linear ODEs associ- 
ated to the classical Painleve equations [3, 4, 5] as well as the similar 
equations associated to the semi-classical orthogonal polynomials have 
matrix dimension 2x2 [8], thus all the relevant spectral curves are (hy- 
per)elliptic. In this section, we discuss the spectral curve for the 3x3 
matrix linear ODE which appears in the theory of the bi-orthogonal 
polynomials [9, 10] with the cubic potentials [19], 




BI-ORTHOGONAL POLYNOMIALS 



8V 



n 



(A) = A n (\)V n (\), * n+1 (A) = R n (\)V n (\), 



n 



u n (\)v 



n 



v n (\)v 



n 



W n (\)V. 



(4.1) 



dx 



dy 



dt 



where 



A n (X) 
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More details can be found in [19] where the Riemann-Hilbert problem 
for \l/ n (A) is formulated. The fixed singularities of the relevant com- 
pletely integrable system correspond to the infinite values of the defor- 
mation parameters x, y, t as well as to t — 0. In the case we are interest- 
ing here, the formal monodromy exponents at A = oo, which is the only 
singular point for the A-equation in (4.1), are equal to n, —n/2, —n/2. 
The asymptotics of \I/ n (A) as n — > oo is of particular importance for 
the theory of coupled random matrices. The scaling changes (2.6) with 
x = — 1/Voo = —1/3 and Remark 2.4 imply 

A = r 1/3 C, x = 5- 2/3 (x + 5 Xl ), 

y = 5- 2 / 3 (yo + 5 yi ), t = S^fa + 5h), n = 5~\ (4.2) 
and yield the system (2.8) with the spectral curve 

F(C, n) = ^ - tl( 3 + fi 2 ( 2 + x /i 2 + yotfc 2 - (f jj - 1R- 

- fiD, - (D 2 -D 3 + 0{5) = 0. (4.3) 

Generically, this curve has 10 first order branch points and therefore, 
via the Riemann-Hurwitz formula, has genus g — 3. Because the mon- 
odromy data for \l/ n (A) are independent from n,x,y,t, the curve (4.3) 
for a generic solution of (4.1) satisfies (3.5), see Remark 3.3. By Theo- 
rem 3.6, given xo,yo,to, system (3.5) uniquely determines the modular 
parameters Dj, j = 1, 2, 3. 

The analysis of (3.5) is significantly more involved then the similar 
analysis of the elliptic curves associated to the classical Painleve equa- 
tions, see [18]. Here, we present the results in the numeric study of the 
integral sub-manifold for (3.5) parameterized by x G C as y = x and 
t = 1 based on the use of MATLAB 6.1 package. 

The graph on the complex ^o-plane shown in Figure 1 separates 
the regions with different topological properties of the relevant Stokes 
graphs. Namely, some of the cycles £j existing in the neighboring re- 
gions collapse at the points of their common boundary. Our numeric 
study suggests that, at the points of the very central triangular domain 
in Figure 1, the curve (4.3) subject to (3.5) has genus g = 0, and the 
relevant Stokes graph is consistent with the Riemann-Hilbert problem 
data of [19]. The latter observation implies the applicability of (3.5) 
to the asymptotic study of the ^-function for the bi-orthogonal poly- 
nomials. In particular, the typical configuration of the branch points 
implied by (3.5) suggests that, for xq ^ 0, the asymptotics of \I/ n (A) 
involves, in certain domains of the complex A-plane, the exponential, 
Airy and parabolic cylinder functions. For x at the boundary of the 
central triangular domain in Figure 1, the asymptotics involves also the 
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Figure 1. The projection of the integral sub-manifold 
yo = x Q , t = 1 for (3.5) on the x -plane 

\I/-function associated to the Painleve first transcendent. For x = 0, 
besides exponential and Airy functions, the asymptotic description re- 
quires also a third order special function. For the values of xo beyond 
this triangular domain, the relevant Stokes graphs seem not consistent 
with the Riemann-Hilbert problem data of [19]. Therefore it is unlikely 
that, for x beyond the central triangular domain, the system (3.5) can 
be applied to the asymptotics of the bi-orthogonal polynomials. The 
detailed description of the asymptotics of the bi-orthogonal polyno- 
mials, however, is out of the scope of the present paper and will be 
published later elsewhere. 
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